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Abstract—The paper presents an exoatmospheric near-optimal explicit guidance algorithm for a satellite 
launch vehicle following a two- or three-dimensional trajectory. The unified steering law has been 
developed in a vector form using an equivalent gravitational field in an inertial coordinate frame. The 
resulting steering law follows a linear tangent law in a canted plane for a three-dimensional trajectory. 
The principle is based on the concepts of optimal transfer between non-coplanar orbits. The guidance 
parameters are determined by solving three simultaneous algebraic equations envolving thrust integrals 
which are computed recursively. Two constant gravity related vectors are defined which account for 
change in velocity and position due to the spherical Earth. The solution is obtained by a differential 
correction method after finding the required partial derivatives analytically. A predictor corrector 
approach is suggested where the trajectory is predicted using the Encke’s method. This also enables 
consideration of other effects such as oblateness. The simulation results for typical multistage launch 
vehicles show that the guidance algorithm is extremely accurate, robust and yet simple enough for on-board 
implementation. It can be used for a variety of missions. 


1. INTRODUCTION 

Closed-loop guidance is essential in upper stages of a 
launch vehicle to reach a specified orbit with min¬ 
imum error. Many schemes have been suggested and 
used during the past three decadesfl]. Most of the 
earlier schemes which are still in use were tailored to 
work under severe hardware limitations of on-board 
processors. With the present level of miniaturization 
and computational power, one can aim at evolving a 
more elegant, optimal and reliable closed-loop guid¬ 
ance scheme which gives higher accuracy, allows 
flexibility of use in a variety of missions, and permits 
fast computation of guidance parameters. This has 
been the objective of the present study. 

The exoatmospheric guidance schemes can be 
classified under two major groups: velocity-correlated 
guidance, and explicit guidance schemes. The 
velocity-correlated guidance schemes[2] take the help 
of the nominal trajectory for a given mission for 
finding the guidance coefficients. Considerable pre¬ 
processing is involved in such schemes and the 
flexibility is limited. Explicit guidance methods are 
based essentially on solution of two-point boundary 
value problems. These schemes are being used widely 
in different forms. In most of the explicit guidance 
schemes available today, steering laws in pitch and 
yaw are derived separately after decoupling the equa¬ 
tions in pitch and yaw planes[3-6]. This requires the 
simplifying assumption of yaw angle being quite 
small. Such decoupling leads to inaccuracy and loss 
of optimality particularly when the pitch and yaw 

tPaper IAF-87-359 presented at the 38th Congress of the 
International Astronautica! Federation, Brighton, U.K., 
10-17 October 1987. 


angles are not small. Use of numerical, rigorous 
calculus of variation methods[7-8] also has not been 
popular because of excessive computation load on 
on-board processor and a possibility of divergence. A 
simplified model is solved[9] semi-analytically with 
simultaneous numerical integration of thrust and 
gravity. However, this approach does not reduce the 
complexity of the calculus of variation problem 
significantly. In Ref.[10j the problem of high yaw 
manoeuvre is addressed and a bilinear tangent law 
for yaw and a linear tangent law for pitch are 
suggested separately. 

This paper presents a near optimal explicit guid¬ 
ance scheme for launch vehicles which follow three- 
dimensional trajectories. The presence of landmass 
constraints and the location of the launch station 
along with the need for vehicle performance opti¬ 
mization sometimes necessitate the use of such 
trajectories, which demand large pitch and yaw 
manoeuvres. A linear tangent guidance law derived 
here in vector form for such missions overcomes the 
shortcomings of the existing approaches. 

The approach adapts some concepts[l 1] of optimal 
transfer between orbits under the assumption of 
uniform gravitational field to find the plane of cor¬ 
rection. The scheme also accounts for the effects of 
actual gravitational field but retains the simplicity of 
the flat-Earth (uniform gravity) approximation by 
defining two gravity related constant vectors corre¬ 
sponding to changes in position and velocity due to 
gravity. The effect of thrust is estimated in terms of 
changes in velocity and position. The parameters 
used in guidance logic, which influence the thrust 
effect, are modified so as to result in a trajectory 
allowing a satellite to be injected into the desired 
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orbit. For predicting the trajectory, Encke’s method 
from orbital mechanics is adapted. For faster com¬ 
putation, the reference trajectory is assumed to obey 
Keplerian relations and is computed analytically. 
This approach facilitates inclusion of oblateness 
effect and other small forces, resulting in a very 
accurate trajectory and helps in finding a near 
optimal solution for vehicle guidance. 

The simulation results show that the algorithm is 
robust and can be used for launch with very high 
accuracy for a variety of missions. 

2. DEVELOPMENT OF ALGORITHM 

Development of the explicit closed-loop guidance 
algorithm requires the determination of the steering 
law in a parametric form. For this purpose first the 
effects of thrust and gravity are estimated separately 
in terms of changes in velocity and position. Then 
three guidance parameters and two gravity related 
vectors are determined which completely define the 
thrust steering programme and the thrust cut-off time 
to meet the desired injection conditions. 

2.1. Coordinate systems 

For the formulation of the algorithm, three right- 
handed rectangular coordinate systems have been 
used; 

(i) The Launch Plane Inertial (LPI) Geocentric 
Coordinate frame OX, Y, Z, shown in Fig. 1(a) is the 
fundamental frame in which all vectors are expressed. 
It has OZ, pointing outward along the launch point 
frozen at take-off and OX, parallel to the launch 
azimuth. 

(ii) Geocentric Orbital Coordinate frame OX G Y 0 
Z Q shown in Fig. 1(b) is defined with OX D pointing 
outward along the present position vector projected 
on the orbital plane, and OY„ perpendicular to the 
desired orbital plane opposite to the angular momen¬ 
tum vector. The unit vector in OY 0 direction t v0 can 
be computed in terms of orbital elements as follows: 

I ro = [4 ] [ - sin (' sin Q, sin i cos fl, — cos (] T , (1) 

where [/(] transforms f ro to LPI frame from an 
equatorial one. 

(iii) In correction Plane Coordinate frame 
OX c Y c Z c shown in Fig. 2, OX c and OY c form the 
plane which contains r thR OY c is along — P r , a costate 
vector related to change in position. 

2.2. Optimal steering law 

In inertial coordinate frame, the basic equations of 
motion of the launch vehicle for its flight in vacuum 
are: 

v = F/m u + g (2) 

i = v (3) 

where v, r are the velocity and position vectors; Fand 
m are instantaneous thrust and mass of the vehicle; 




MERIDIAN 
AT LAUNCH 

(b) 

Fig. 1. (a) Launch Plane Inertial Coordinate frame, (b) 
Geocentric Orbit Coordinate frame. 

g is the gravitational acceleration and u represents the 
unit thrust vector. 

The optimization criterion chosen here is the min¬ 
imization of time of burn for the rocket (or the last 
stage in the case of a multistage rocket) where 
thrust-time profile(s) (and the burning sequence of 
the different stages) is (are) specified, the problem 

LOCUS OF POINT P 
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of minimum time transfer under the assumption of 
uniform gravity is well documented[l 1]. This can be 
used with appropriate modification to develop the 
guidance logic. The solution can be written in a 
vector form as: 


Pv= -Pr. 

P r = 0. 

Thus, 

Pv(0=-Prf+C, (4) 

where p v is the Lawden's prime vector which is a 
costate vector associated with velocity and p r is a 
constant vector defining the costate associated with 
position. Another constant vector is c. Equation (4) 
shows that p v (r) is contained in a plane generated by 
c and p r , as shown in Fig. 2, where OA represents c. 
After resolving c into two components, along p r and 
perpendicular to it, let the vector perpendicular to P r 
be p, along OB with the associated unit vector f„. Let 
AB = — p r r c where r c is a constant having unit of 
time. Hence, 

C 4“ Prt c . 

On using this equation in eqn (4), one gets, 

Pv(0/<c = J re + Po(t ~ * c ), ( 5 ) 

where f,, c is the unit vector in the direction of —p,, 
and pa = |p r |/r c . The unit vector u in the direction of 
thrust can be obtained from eqn (5) by dividing the 
right-hand side by its magnitude, noting that T vc and 
i K are orthogonal. 

U(/) = [L + lycPtO - T c )]/[1 +p 2 o(t ~ t c ) 2 ] ,/2 . (6) 

This constitutes the optimum steering law. It can be 
seen from eqn (6) that the angle V, which the thrust 
vector makes with I IC can be expressed as 

tan F = p„(t -t c ). (7) 

Hence, the tangent of the optimal steering angle 
follows a simple linear tangent law under uniform 
gravitational field. 

2.3. Terminal constraints 

The constraint relations at injection can be com¬ 
puted by first computing the desired injection state 
and then equating it with the integrated trajectory 
following a steering law up to thrust cut-off time. 

For orbital injection, generally five out of six 
possible parameters are specified. They are R D , V D , 
y D and the desired orbital plane defined by unit vector 
normal to the orbit. The injection state then can be 
defined in terms of range angle 0 X and unit vectors, 
! to and (defined in Section 2.1) as: 

V T = Vb&o sin(y D - Or) + L 0 cos(y D - 0 T )], (8) 

f t = 7?d ffxo cos(0 x ) + I zo sin(0 T )]. (9) 

Alternatively, the injection constraints can be defined 
in terms of orbital elements, namely a, e, 6 , i, fi, 


where i and C2 define the orbit normal, to is uncon¬ 
strained. The corresponding values of R D , V D , and y D 
can be determined using the following relations: 


R d = u(l — e 2 )/(l -1- e cos 0), 

(10) 

V D = [p(2/R B -\/a)Yf 

(11) 

? D = cos _1 [{/ta(l — e 2 )Y l2 l(R D V D )}. 

(12) 


This required state at injection leads to an equality 
constraint obtained by integrating the equations of 
motion up to T go using u(/), which in turn is defined 
in terms of p$, z c , t xc and I >1C . It may be noted that 
there are only five constraints specified, hence injec¬ 
tion can take place anywhere in the orbit. As the 
steering law is derived under the assumption of 
uniform gravitational field, it is valid for zero gravity 
field also, with appropriate change in the end con¬ 
straints. The transformed end conditions allow one of 
the position components, along which total velocity 
constraint is met, to be unconstrained. It can be 
shown[12] that this corresponds to the direction of I xc 
(Fig. 2). Since the velocity constraint is met in the 
direction of I w , the net change in the velocity over the 
entire path in I re direction should be zero. This 
relation is used to define r c . Since I vc and f re are 
orthogonal, one is left with five independent param¬ 
eters associated with u(r), (three defining ! IC and ! K , 
one each for p,■, and t c ). 


2.4. Effects of thrust 

The changes in velocity and position due to thrust 
are estimated along with the corresponding changes 
due to gravity with assumed values or first approxi¬ 
mation of the steering parameters. The values of the 
parameters are then modified in such a way that 
the resulting thrust direction programme guides the 
launch vehicle from the present state to the desired 
injection state. 

The required change in velocity to be imparted by 
thrust is given by 


(13) 


where 


gMO} dr. 


(14) 


The estimated change in velocity imparted by 
thrust for a specified set of steering parameters and 
F go is given by 

rr , 0 


F(t)u(t)/m{t)dt. 


(15) 


On substituting eqn (6) in to eqn (15) one obtains 


r i 

Jo 


F(t)l[m(t){\+p](t-x c m&t 


+ 1 


■i 


F{t)pd{t - t c )/[m(r) 


x{l+^(r-r c ) 2 }n 


(16) 



118 


S. K.. Sinha et al. 


As explained in Section 2.3, the velocity constraint is 
to be met by the component of thrust in I tc direction 
only, therefore, 


V,hA = I vc 

and 
0 = 


Fmrn(t){l+pUt-i c ) 2 }' !2 ]dt (17) 


F(t)Pi(t -x c )/[m(t) 

x{l + pj(t — t c ) 2 }' /2 ] dr. (18) 


Similarly, considering the position constraint, one 
gets the required change in position due to thrust: 

r thR — r T -r 0 -v o r g0 -r g , (19) 

where 


r 0 = 


•r g0 rt 
o Jo 


g{r(j)} ds d t. 


r thA — I vc 


0 J 


F(s)/[m(s) 


x {1 +p^(j -r c ) 2 } l/J ]ds d t 


+ t, 


'T, „ 

Jo . 


F(s)p 6 (s -T c )/[m(s) 


X {1 -t-/>^(s — T c ) 2 } 12 ]dv dr. 


( 20 ) 


The estimated change in position is given by: 
rr, o 


( 21 ) 


For the solution, the following constraints must be 
satisfied. 

V,hR = V,hA, (22) 

f thR = r ihA • (23) 

Equation (18) is solved to determine the value of x c . 
Since closed-form solutions of eqns (17), (18) and (21) 
are not available, some simplifications are made to 
make them integrable. For this the following relation 
is used. 

[I + pj(t — A-) 2 ]" 1/2 = 1 ~2Po(‘ - A) 2 

+ 1 /4(r-t c ) 4 ----. (24) 

After selecting a finite number of terms of the above 
series and using eqn (24) in eqns (17), (18) and (21), 
one gets equations of the form: 


» T 

F(t)t"/m(t) dt, 

Jo 

(25) 

j F(s)s n /m(s)ds d/. 

Jo Jo 

(26) 


These integrals can be solved analytically for con¬ 
stant thrust cases. For this scheme, the integrals are 
solved in a recursive way [12] making the approach 
computationally efficient. Equations (17), (18) and 
(21) can be written in algebraic form in terms of the 
integrals I Sn and I Dn . For example, on using the first 
two terms of expansion from eqn (24), one gets the 


following equations: 

l ’thA = 7so — iPo(Is2 ~ 2r c /si + t~/ so ) (27) 

0 = C4i — A / s „) — \pii(fsi — 3t c /s: 

+ 3f 2 / si — t^/ S o ) (28) 

r thA = [/»„ — 2Pri(^n; — 2r c / D | + r~/ Do )]~ 

+ PiiK^Di — A/d.,) — jpl. >Ud} ~ 3r c /p2 
+ 3T c 2 / DI -r t '/ Do )] : . (29) 

Although two terms of eqn (24) for expansion have 
been considered here, one can take as many number 
of terms as required for the desired accuracy. 


2.5. Gravity effect modelling 

The effect of gravity depends upon the trajectory 
followed by the launch vehicle, which in turn depends 
upon the value of the steering parameters. In con¬ 
formity with the assumptions made in the derivation 
of the steering law, the gravity vector has to be taken 
to be uniform over the entire path. In order to 
represent the effects of the actual gravity on position 
and velocity vectors, two gravity related constant 
vectors are defined below: 


g 

and 

g, 


-n 


g[r(0] d t > 


dt 


g{r(s)} ds dr / 



ho 

ds dt 


0 . 

“I 

o 


On using eqns (14) and (20), one gets, 

gv — v g /7- g0 

and 

g r = 2r g /T" 0 


(30) 


(31) 


(32) 

(33) 


By choosing these two gravity related constant vec¬ 
tors, the simplicity of uniform gravity is maintained 
without sacrificing the accuracy of trajectory gener¬ 
ation. For computing v g and r g , Encke’s approach of 
integrating the perturbations has been adapted to 
meet the fast computation requirement of the 
closed-loop guidance. The reference trajectory used 
corresponds to the Keplerian relations for two-body 
motion. The other forces acting on the vehicle are the 
thrust and a small force due to oblateness of the 
Earth. Encke’s approach requires acceleration due to 
these perturbative forces alone to be integrated. 

Let r E be the position on the reference trajectory 
and r be the position on the actual trajectory which 
the vehicle follows when all forces act on it. The 
equations of motion can be written as: 

r= -pr /r 3 + Ar, hrust + Af ob , atc (34) 

where A¥ thrust and Ar oblale are the accelerations due to 
thrust and oblateness of the Earth respectively. 
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Let 

% — r — r E . (35) 

On using eqns (34) and (35), one gets, 
l = Hl r \[r{\ -(r E /r) 3 } -£] + Af thrust + Ar oblate (36) 
or 

\ = r E + $)f{q) - %\ + Ar lhrust + Ar oblate (37) 
where 

f(q)=\-(\+2qY^ (38) 

and 

q = {^(T E + (,/2)}/rl (39) 

Equation (37) can be written in a concise form as: 

l = Ar Encke + Ar thrust + Ar oblale . (40) 

r E and v E can be related to the state at t a by a 
transition matrix: 


r E (r) 

, v e(0. 


= #(*,*.) 


^*E (t 0 ) 
_ V e('o). 


(41) 


The solution for r E (r) and v E (?) are obtained analyti¬ 
cally. Integration of Ar thrust can the performed using 
Section 2.4. Ar Encke and Ar oblate need to be integrated 
numerically. Since these vectors are of very small 
magnitude as compared to that of r E , large step-size 
can be used for integration without sacrificing the 
desired accuracy. The expression for Ar oblale is: 


Afobime = ~/t[3/2]7 20 (7? c /r ) 2 


X {1 - 5(u.-r) 2 /r 2 }]r/r 3 - 3pJ 20 (RJr) 2 
x (u.T)u./r 3 (42) 


where J 20 is the second zonal harmonic term, u. is the 
unit vector in the direction of the North Pole and 7? e 
is the radius of the Earth. For numerical integration 
of Ar Encke and Ar oblate , Simpson’s scheme is found to 
be satisfactory. 

Integration of eqn (40) along with eqn (41) pro¬ 
vides the prediction at injection. With the above 
formulation, the predicted injection state is generally 
quite close to the desired one. This prediction helps 
in the modification of the guidance parameters to 
meet the injection constraints accurately. It may be 
noted that g v and g r computed from the predicted 
injection state converge to the correct value very fast 
as they represent the time averaged value over the 
entire path. 

The vectors g v and g, are computed as follows: 


gv = (Vr-v 0 -v thA )/r g0 (43) 

gr = 2(r t - r 0 - v r go - r thA )/ T\ 0 (44) 

where vj and rj are the predicted velocity and posi¬ 
tion vectors at the injection using Encke’s method. 


2.6. Guidance parameters 

The parameters I tc , T lt , p „, t c define the steering law 
[eqn (6)] completely. These are calculated after deter¬ 


mining the parameters T go , p„ and 0 r which are found 
by iteration after converting the vector eqns (22) and 
(23) into the corresponding algebraic scalar equations 
and applying a differential correction approach. The 
scalar algebraic equations corresponding to eqns (22) 


and (23) are: 

^1 =v thR -v lhR -{(v thA -I Jtc ) 2 } (45) 

F- = r thR ' r th R — ( r thA ' ) 2 — ( r thA ' Y (46) 

F3 = v lhR • r lhR — v thA • r thA (47) 

^ 4 = v thA I K . (48) 


r c is determined using eqn (48). At the correct value 
°f 7go> Po a °d ^t> ^T, F2, F3 should be simulta¬ 
neously zero. Any small change in their value results 
in an error in F t , F 2 , Fj according to the following 
relations: 

6F\/dT lo dT go + dFl/dpt d P(j 

+ dFl/d6 T dQj = dF\, (49) 

dF2ldT go dT go + dF2ld Pl idp 0 

+ 5F2/86 t d0 T = dF2, (50) 

dF3/dT go dT go + dF3/8 P6 d Pli 

+ dF3/dO T d0 T = dF3, (51) 

where the partial derivative can be found analytically 
by differentiating the right-hand side of eqns (45), 
(46) and (47) with respect to T g0 , p,j and 0 T . From 
these linearly independent algebraic equations, the 
corrections needed in the estimates of T g0 , p,j and 0 T 
to make Fi(i = 1,2,3) zero can be found from the 
errors in Fi by matrix inversion. The scheme exhibits 
good convergence. Generally, it takes only one 
or two iterations for the solution except at the 
start of guidance where it may take up to 5-6 
iterations depending upon errors in initial estimates 
of guidance parameters. The convergence of T go 
[i.e.|d7' g0 /7’ g0 | < £ (a small number)] ensures the 
convergence of p,, and 0 T as well. 

With these, the value of v lhR and r, hR can be found 
by using eqns (13) and (19) respectively. I vc and I, c are 
then determined as follows: 

L = unit (v thR ) (52) 

L, = unit x r lhR ) x I tc . (53) 

Thus, all parameters in eqn (6) defining the steering 
law are fully determined. 

2.7. Terminal processing 

In the vicinity of the injection point, most of the 
guidance schemes tend to instability. This originates 
from the terminal constraint which requires zero 
injection error, and which may demand infinite 
control gains and very fast turn over rates close 
to injection. Therefore, one generally switches over to 
another guidance mode. 

In the present case, the position constraint imposed 
by eqns (46) and (47) are relaxed once T g0 is less than 
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r min , a predefined time based on simulation. There¬ 
after, only one equation, namely eqn (45) is solved for 
T go , keeping the injection state and p„ the same as 
determined earlier. 


2.8. Multistage rocket 

For a multistage rocket the guidance algorithm 
needs some modification. The computation of v lhA 
and r thA are different, as one has to compute the total 
changes in velocity and position imparted by the 
current and remaining stages. Also, while solving the 
double integrals, one should include the con¬ 
tributions by the corresponding single integrals at the 
end of the previous stages. Other computations are 
identical to the ones corresponding to the single 
stage. 

As it is simpler to compute the thrust integrals of 
eqns (25) and (26) between the limits of zero and a 
definite time, the value of t c has been correspondingly 
modified for the remaining stges. Thus 

A f+( =*c, ~ At nj < 54 > 

where r c . is z c for the current stage and A t l+J is the 
time difference between ignition time of (/ + j) stage 
and the current time. This way t c for all the remaining 
stages can be defined. 

The expression for v thA can be written as, 


T ‘ FAm, 
o 

x{l+^(/-T c ,) 2 }' ,2 ]dr (55) 

where T, is the remaining burn time for the ith 
stage. t c for the current stage is computed using the 
following relation: 



0= I 

/ = current stage J 


FiPoit - F)/K 

x (1 + Pt(t — T c .) 2 } 1/2 ]dr (56) 


t c for all stages are expressed in terms of t c for the 
current stage, which is then found using eqn (56). 

The expression for r lhA can also be written in a 
similar manner on using the following relation: 



i = current stage 



m t 

F i (s)/[m(s) 

0 


X {1 +pj(s — t C( ) 2 } 1?2 ] ds dr 


+ >„ 


FMUmU^l+pUt-xJ 2 }' 12 ]^ 


+ 1, I 

/ = current stage 

X {1 +/^(i -T Ci ) 2 } ,;2 ]di dt 
% T< 

+ t<, ' 


F,(s)pg(s -T c .)/[m(s) 


F,m>n{t){\ + p](t -x Cl f} xll ]dt 


( 57 ) 


where r r is the remaining time for the thrust cut-off 
after the end of the /th stage. 

For computation of gravity and prediction of state 
at the expected injection, Encke’s method is used. The 
trajectory is corrected at the end of each stage using 
the v lhA and r lhA corresponding to the particular stage 
as given in Section 2.5. 

The guidance parameters are determined in a man¬ 
ner identical to the one given in Section 2.6. 

3. SIMULATION RESULTS 

Performance studies are carried out for a represen¬ 
tative multi-stage launch vehicle which employs a 
3-dimensional trajectory due to land-mass constraint, 
for launching a satellite into a sun-synchronous orbit 
of altitude 900 km. The closed loop guidance is 
operational for nearly 1000 s including a long coast 
before the final stage ignition. In total, it covers a 
range angle of nearly 40°, requiring nearly 120 of 
manoeuvre in pitch and nearly 20 c manoeuvre in 
yaw. The update of guidance parameter is done for 
every 5 s and the steering angle is computed based on 
determined guidance parameters at every 0.5 s. The 
position constraint is relaxed if T i0 <5s. For esti¬ 
mating the thrust effect, the first two terms of series 
expansion are retained. 

The salient features of the results are shown in 
Table 1. The results for a nominal flight shows very 
good accuracy at injection. The error in semi-major 
axis is clOcm and that in eccentricity is around 
10 7 . The position error and errors in apogee and 
perigee are all less than 1 m. Figure 3 shows the pitch 
and yaw angles for the duration where guidance was 
operational. During the coast, steering is ineffective, 
hence, these angles have not been shown. 

This high accuracy has been possible due to accu¬ 
rate prediction of the trajectory. The prediction re¬ 
sults are shown in Table 2. It shows that even with 
no update of guidance parameters for the last 150 s, 
the error in semi-major axis is <16m and error is 
eccentricity is around 10 ~ 5 . The prediction error has 
been found by first computing the injection state by 
numerical integration using the vehicle's position and 
velocity at some previous time with the steering 
parameters computed at that time and then com¬ 
paring the injection condition with the desired one. 

Many cases of perturbations have been simulated 
to test the robustness of the algorithm. Typical results 
of the following cases have been shown in Table 1. 

Case 1. This case has been simulated to test the 
explicitness of the algorithm by varying the initial 
position and velocity at the initiation of the guidance. 
For this, each component of the position vector is 
changed by 50 km and each component of velocity 
vector is changed by 50 m/s. This does not affect 
adversely the initial convergence or the final accuracy 
of injection. The change in burn-time is 17.8 s. 

Case 2. While considering the variation in the 
system parameters, the most dominant effect of per- 
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Table 1. Simulation of injection error—Sun-synchronous orbit 


Case 

No. 

Description 

Change in 
burning-time 

(s) 

Error in 
eccentricity 

Abs. error 
in semi¬ 
major axis 
(m) 

Abs. error 
in orbit 
normal 
(deg) 

Abs. error 
in inj. 
position 
(m) 

Abs. error 
in inj. vel. 
(m/s) 


Nominal performance with no perturba¬ 
tion 

— 

1.24 x 10‘ 7 

0.086 

7.54 x 10“ 6 

0.955 

2.88 x 10~ 7 

i. 

Change in initial positions: 50 km in each 
components Change in initial velocity: 
50 m/s in each component 

17.8 

1.05 x 10"‘ 

0.10 

6.0 x 10- 5 

7.63 

1.61 x 10" 5 

2. 

Change in specific impulse: penultimate 
stage, 1%; final stage, 0.75% 

6.36 

1.74 x 10 5 

0.33 

9.9 x 10 4 

4.39 

0.133 

3. 

Change in specific impulse for all stages, 
5% 

57.12 

1.33 x 10 4 

116.2 

7.54 x 10-’ 

297.04 

0.675 

4. 

Change in structure weight: penultimate 
stage, 3%; final stage, 1.5% 

8.94 

3.6 x 10 s 

10.05 

2.05 x 10- 3 

266.17 

5.52 x 10- 3 

5. 

All perturbations included: Change in 
specific impulse: same as case 4 change 
in propellant wt: penultimate stage, 0.5% 

19.1 

7.86 x 10 s 

52.71 

4.48 x 10 3 

577.24 

1.63 x 10“ 2 



Fig. 3. Pitch and yaw angles for a Sun-synchronous mission. 


turbation is due to variation of specific impulse of the 
propellant used in the motor. This case is simulated 
to test the robustness of the algorithm under 3 a 
variation of the specific impulse, which is taken as 1 % 
in the penultimate stage and 0.75% in the final stage. 
The error is semi-major axis is 0.33 m and error in 
eccentricity is 1.7 x 10 5 . The change in burn-time is 
small. 

Case 3. This case also has been simulated with a 
reduction in specific impulse in all the stages by 5% 
as compared to the nominal one. This high value of 
specific impulse variation has been considered to 
highlight that the algorithm is very robust. The 
convergence takes nearly 3 iterations most of the 
time. The error in semi-major axis is around 120 m 
and eccentricity error is 1.3 x 10 \ which is well 
within the specification of the injection accuracy of 
such missions. The penalty is significant increase in 
the burn-time of the last stage. 

Case 4. The source of perturbation other than the 
variation in the mass of the vehicle. This has been 
simulated by considering 1.5% of higher structural 
weight in the final stage and 3% of higher structural 
weight in the penultimate stage. The resulting error 
in semi-major axis is around 10 m, and error in 
eccentricity is 3.6 x 10 -5 . This result shows that the 
algorithm is quite robust. 

Case 5. The robustness of the algorithm has been 
further tested under the extreme perturbation case 


Table 2. Prediction accuracy 


Parameter for which 
error is computed 


Prediction duration 


150 s 

250 s 

400s 

Semi-major axis (m) 

15.7 

-63.3 

-602.8 

Eccentricity 

2.27 x 10 s 

3.79 x 10“ 5 

1.53 x lO" 4 

Position X (m) 

-23.27 

-41.37 

-9.26 

Position Y (m) 

125.57 

214.40 

456.95 

Position Z (m) 

103.69 

172.63 

582.96 

Velocity X (m/s) 

-2.53 x 10 3 

-2.89 x 10~ 3 

-2.60 x 10~ 2 

Velocity Y (m/s) 

1.59 x 10“ 3 

3.30 x 10“ 2 

5.33 x 10“' 

Velocity 2, (m/s) 

-6.96 x 10" 3 

2.29 x 10“ 2 

2.54 x 10 ' 
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Table 3. Case study of geostationary transfer orbit injection 


Error in perigee height 

13.01 m 

Error in apogee height 

-17.15m 

Error in semi-major axis 

-2.07 m 

Error in eccentricity 

6.4 x 10“’ 

Error in injection position: 

Component X 

-0.89 m 

Component Y 

0 m 

Component Z 

0.16 m 

Error in injection velocity: 

Component X 

—1.85 x 10~ 6 m/s 

Component Y 

0 m/s 

Component Z 

1.35 x 10 6 m/s 


where perturbations due to variation in specific im¬ 
pulse, structure weight and propellant weight all have 
been considered simultaneously as follows: increase 
in structure weight: 3% in the penultimate stage, 
1.5% in the final stage; decrease in specific impulse: 
1% in the penultimate stage, 0.75% in the final stage; 
decrease in propellant weight: 0.5% in the pen¬ 
ultimate stage. The final error in semi-major axis is 
round 50 m and error in eccentricity is around 
8 x 10~ 5 with 19.1 s increase in burn-time. 

All these results prove that the algorithm is robust 
under various perturbations possible during flight 
and is capable of injecting the satellite into the desired 
orbit very accurately. 

Performance of the scheme has also been tested for 
a geostationary transfer orbit mission which employs 
a planar launch trajectory. The launch vehicle for 
which simulation studies have been made has three 
stages and closed-loop guidance is operational in the 
last two stages. The result in Table 3 shows that the 
accuracy of the injection is good. The error in 
semi-major axis is only 2 m and error in velocity is 



less than 2 x 10 _6 m/s. The position error at injection 
is less than 1 m. Figure 4 shows the pitch angle 
and pitch rate during the period when guidance is 
operative. 

4. CONCLUSION 

An accurate, efficient and truly explicit closed-loop 
guidance algorithm has been formulated in a vector 
form using Launch Plane Inertial Coordinate frame 
which works well for a launch vehicle following 2- 
and 3-dimensional trajectory even with relatively 
large range angle and long operation time. 

The algorithm is quite flexible in its computational 
load demand. The number of terms considered for 
thrust effect estimation can be properly chosen to suit 
the on-board computational power and a trade-off 
can be made between prediction accuracy and com¬ 
putational load. The simultaneous computation of 
guidance parameters avoids the need of separate 
modules to compute time-to-go, range angle and rate 
of turn of thrust angle and ensures fast and true 
convergence of guidance parameters. The accuracy of 
predicted trajectory using Encke's method allows the 
guidance parameter update cycle to be longer and at 
the same time helps in the getting a steering solution 
which results in an optimal and accurate guidance 
scheme. 

The inertial navigation coordinate frame used in 
the formulation avoids all coordinate transformation 
requiring computations of trigonometrical functions. 

With its high accuracy and flexibility, the algorithm 
can be used for a variety of missions. Here its 
applicability in the case of a 3-dimensional launch to 
reach a Sun-synchronous orbit and a 2-dimensional 
launch to reach a transfer orbit for geostationary 
missions in demonstrated. 
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APPENDIX 

Nomenclature 

a, e, i, co,Cl,6 = semi-major axis, eccentricity, inclina¬ 
tion, argument of perigee, longitude of 
the ascending node and true anomaly 
of the resulting orbit at injection 
F,m= instantaneous thrust and mass of the 
vehicle 


(9 = x, y, z) = unit vectors in correction and orbital 
coordinate frames 

Pi = a guidance parameter corresponding 
to rate of change of thrust angle 
r g ,v g = change in position and velocity due 
to gravity computed from the present 
time until injection 

r o> v o - present position and velocity obtained 
from navigation system 
r T ,v T = desired position and velocity at 
injection 

r thA • r thR = change in position estimated and 
required to be imparted by thrust 
R D , ^d> Vd ~ desired distance, velocity and flight 
path angle at injection 
T g0 = time-to-go, the thrust cut-off time 
measured from the present time 
T'min = the value of T m after which position 
constraint is relaxed 

v ihA> v thR = change in velocity estimated and 
required to be imparted by thrust 
6j = range angle in the orbital plane 
measured between the projection of the 
present position vector and estimated 
injection location 

r c = time parameter in the steering law 
p = Earth’s gravitational constraint. 
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